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Abstract 



We review the theory of polarized deep inelastic scattering in light of the most recent experimental 
data. We discuss the nucleon spin decomposition and the Bjorken sum rule. The latter is used for 
extraction of a s {M%) = 0.116±g;^ and as a test case for a new method of analyzing divergent 
perturbation series in QCD. 



Analysis of Polarized Structure Functions 

Formalism 



The basis for our discussion will be the two spin- 
dependent structure functions G± and G 2 '- 
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In the parton model, these structure functions scale 
as follows in the Bjorken limit x — Q 2 /2M^v fixed, 

* Based on plenary talks presented by M. Karliner at the 
Workshop on Deep Inelastic scattering and QCD, Paris, April 
1995 and at the Workshop on the Prospects of Spin Physics at 
HERA, DESY-Zeuthen, August 28-31, 1995. 



We will discuss the scaling structure function g 2 later 
on, focussing for now on g\ , which is related to the po- 
larized quark distributions by 



2^2 e «[9T0) 
1 



Ql(x) +<7t0) - 



for comparison, the unpolarized structure function F2 is 
given by 

F 2 {x) = ^2e 2 q x[q^{x) + q l {x)+q 1 {x)-q l {x)] (4) 

so that the polarization asymmetry A\ may be written 

as 

cr i/2 - C3/2 



At 



(5) 



Cl/2 + °3/2 

in the Bjorken limit, where CT1/2 and 03/2 are the virtual 
photon absorption cross sections. We will discuss later 
the Q 2 dependences of the above formulae, as well as 
the transverse polarization asymmetry. 

Much of the interest in the polarized structure 
function 171 is due to its relation to axial current matrix 
elements: 

(p\Al\p) = (p|«7M7s«b) = (pWrJ/mQr - Ql7hQl\p) = 



Aq-S^p) 



(G) 



where q L R = 1/2(1 =p 75) q, is the nucleon spin four- 
vector, and 



Aq = dx[q^ (x) - q x (x) + q } (x) - q l (x)] (7) 
Jo 

Of particular interest is the matrix element of the singlet 
axial current 



(p\Al\p)= A 1-S,(P) 

q—u,d,s q—u,d,s 

(8) 

which is related in the parton model to the sum of the 
light quark contributions to the proton spin. Prior to 
the series of measurements of polarized deep inelastic 
lepton nuclear scattering, information was available 
from charged current weak interactions on some axial 
current matrix elements. For example, neutron beta 
decay and strong isospin symmetry tell us that jlj 



Au - Ad = F + D = 1.2573 ± 0.0028 



(9) 



and hyperon beta decays and flavour SU(3) symmetry 
tell us that |) 



Am - 2As 



(IK 



3F-D 



V3 



a/3 a/3 



0.34 ±0.02 (10) 



(for a recent discussion of the applicability of SU(3) 
symmetry see Ref . |jj and references therein) . Equations 
(S) and ( |l0| ) give us two equations for the three 
unknowns Ait, Ad and As. In principle, a third piece 
of information was available j|, |5| in 1987 from neutral 
current weak interactions. Measurements of vp and Dp 
elastic scattering M indicated that 



As 



-0.15 ±0.09 



(11) 



but this information was not generally appreciated 
before the advent of the EMC data discussed below. 
At present there is a new neutrino experiment under 
way at Los- Alamos which is expected to significantly 
improve the precision of (|ll|) (see Ref. || for a recent 
in-depth analysis). 

In the naive parton model, the integrals of the 171 
structure functions for the proton and neutron 



r?(Q 2 )^ fdxgl{x,Q 2 ) 
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dxgUx,Q 2 ) 



are related to combinations of the Aq. 



T'P = -( -Au + -Ad + -As 
1 2 1 9 9 9 
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The difference between the proton and neutron 
integrals yields the celebrated Bjorken sum r ule § 

r?(Q 2 ) - r?(Q 2 ) = i(A«- Ad) x (1 - a, (Q 2 )M + ... 

(14) 

It is not possible to derive individual sum rules 
for r?'" without supplementary assumptions. The 
assumption made by Ellis and Jaffe in 1973 |[o| was that 
As = 0, on the grounds that very possibly there were a 
negligible number of strange quarks in the nucleon wave 
function, and if there were, surely they would not be 
polarized. With this assumption, it was estimated that 

r 1 1 

/ dxgP(x, Q 2 )= — (4Au + Ad) (1 - a s /TT + ...) = 



(15) 



= 0.17 ±0.01 



It should be clear that this was never a rigorous 
prediction, and was only intended as a qualitative 
indication to experimentalists of what they might find 
when they started to do polarized electron proton 
scattering experiments. 

Perturbative QCD corrections to the above relations 
have been calculated lO-IO: 



^(^Q 2 )-5r^Q 2 )] = ^l.9A|/(x) : 
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f(x) = 1 - x - 3.58x 2 - 20.22x 3 



(16) 
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and 



9 P i in \x,Q 2 ) = 



h{x) = l-x- 1.096a; 2 - . . 



-A£(Q 2 )/i(x) 



(17) 



where x = a s (Q 2 )/ir, and the dots represent 
uncalculated higher orders of perturbation theory, to 
which must be added higher-twist corrections which 
we will discuss later. The coefficients in (|l^),(^7|) are 
for Nf=3, as relevant for the Q 2 range of current 
experiments. 

With these corrections, the Bjorken sum rule is a 
fundamental prediction of QCD which can be used, for 
example, to estimate a value for a s (Q 2 ). On the other 
hand, the individual proton and neutron integrals can 
be used to extract a value of As. 



1.2. The Helen of spin 

Early data on polarized electron-proton scattering from 
SLAC-Yale experiments JlJ, [D| |l6[| were compatible 
with the prediction of equation ( |15| ) within large errors. 
Over a 1000 theoretical and experimental papers were 
launched by the 1987 EMC result 

/ g^(x,Q 2 ) = 0.126 ±0.010 (syst.) ±0.015 (stat.) , 
Jo 

at (Q 2 ) = 10.7 GeV 2 (18) 

which was in prima facie disagreement with the 
dynamical assumption that As = 0. It is worth pointing 
out that the small- x behaviour of g\{x) was crucial 
to this conclusion. The earlier SLAC-Yale data had 
large extrapolation errors, and the EMC data indicated 
behaviour different from that in simple dynamical 
models. They were, however, consistent Q with the 
naive Regge expectation Eq] 



01 0*0 - ^CiX Ql 



(0) 



(19) 



were the «i(0) are the intercepts of axial vector Regge 
trajectories which are expected to lie between and -0.5. 
A fit to the EMC data gave B 
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S = -0.07i°;g for x < 0.2 (20) 



for x < 0.2. 

Using equations (||) , ([l0|) , (|lS|) and the leading order 
perturbative QCD corrections in equation ([T^) it was 
estimated Q that 



Am = 0.78 ±0.06 
Ad = -0.47 ±0.06 
As = -0.19 ±0.06 



(21) 



Strikingly, these determinations corresponded to a total 
contribution of quarks to the proton spin 



AS = Au + Ad + As = 0.12 ± 0.17 



(22) 



which was compatible with 0. This has sometimes been 
called the "proton spin crisis" , but we think this is an 
over-reaction. The result equation (22) was certainly 
a surprise for naive models of non- perturbative QCD, 
but it was not in conflict with perturbative QCD. 
Moreover, shortly after the first data became available 
it was shown |l^] that AS = occurs naturally in 
the Skyrme model, which is believed to reproduce 
the essential features of QCD in the large-iV c limit. 
Alternatively, it was suggested that the U(l) axial 
anomaly and polarized glue might provide an alternative 
interpretation pTj|-|p2[, or a significant suppression of 
the QCD topological susceptibility p3|-[p5| might play 
a key role, which would modify the naive quark model 
predictions. 

In 1972 Richard Feynman wrote ". . . its [the Bjorken 
sum rule's] verification, or failure, would have a most 
decisive effect on the direction of future high- energy 
physics" . On the other hand, we think that the 
verification, or failure, of equation ( p2| ) has only an 
indecisive effect, though a very interesting one. 

1.3. Evaluation of integrals 

Before discussing the interpretation of more recent data 
on polarized structure functions, we first review a few 
points that arise in the evaluation of the integrals F^'™. 
It should not be forgotten that the QCD versions of the 
sum rules are formulated at fixed Q 2 . A generic deep- 
inelastic sum rule in QCD reads 



r(Q 2 ) = r c 



l 



n>l 



a s (Q 2 ) 



m>l V ^ ' 



(23) 

where Y ^ is the asymptotic value of the sum rule for 
Q 2 — > co, the c„ are the coefficients of the perturbative 
corrections, and the d m are coefficients of the so-called 
mass and higher-twist corrections. On the other hand, 
the data are normally obtained at values of Q 2 that 
increase on the average with x Bj as seen in Fig. [IJ. 

It is therefore necessary for each individual experi- 
ment to interpolate and extrapolate to some fixed mean 
value of Q 2 , as indicated by the dashed horizonal line 
in Fig. |l|. The quantity measured directly is the polar- 
ization asymmetry ([s]), which seems experimentally to 



3 




behaviour as x — > 0, so that 



0.0 



0.5 
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Figure 1. In any given polarized lepton-nucleon scattering 
experiment, the range of Q 2 probed is different in different bins 
of the Bjorken variable Xb ■ ■ 



have only small dependence on Q 2 . Therefore experi- 
ments often assume that A\ is a function of x only, and 
then estimate 



9i(x, Q 2 ) 



A 1 (x,Q 2 )F 2 (x,Q 2 ) A 1 (x)F 2 (x,Q 2 ) 



2x[l + R(x,Q 2 )] 2x[l + R(x,Q 2 )} 



(24) 



where F 2 (x, Q 2 ) and R(x, Q 2 ) (the ratio of longitudinal 
to transverse virtual photon cross-sections) are taken 
from parametrizations of unpolarized scattering data. 
Note that these induce a Q 2 dependence in g\ even if 
A\ is independent of Q 2 . 

The possible reliability of the assumption that A\ is 
independent of Q 2 can be explored using leading-order 
perturbative QCD models for Qi(x, Q 2 )- Several such 
studies have been made p6|-|pl|, and they indicate a 
small Q 2 dependence in Ai(x,Q 2 ). This has in turn a 
small effect on the extracted values of gi(x, Q 2 ), which 
is much smaller than the statistical errors of the EMC 
and SMC experiments and is not very significant for 
the E142 and E143 experiments, and is in any case 
considerably smaller than the systematic errors, so that 
it does not yet contribute an important error to the 
evaluations of the 1^* . However, it could become an 
important effect in the future, and both theorists and 
experimentalists should keep their eyes open. 

The old-fashioned assumption of Regge behaviour 
at low x also needs to be checked carefully. The 
leading-order perturbative QCD evolution equations 
for the non-singlet part of the helicity distributions, 
AqNs( x iQ 2 ), l ea d us to expect |52], |33| singular 



1 



Aq NS (x,Q ) ~ C NS exp(A NS a + B N s \np--\n a) 

P 2 

(25) 

where Am a. Bmq and Cnq are some constants and 



x t 




(26) 

and we might expect by analogy with the unpolarized 
structure functions that xq ~ 0.1, Qq ~ 1 GeV and the 
leading-order QCD scale parameter A ~ 0.25 GeV, with 



Ans — 



4^2 



V/33-27V/ 



Bns — 



33 - 27V 



f 



(27) 



where Nf = 3 in the Q 2 range of current experimental 
interest (see also Ref. |34|]). In principle Eq. (25) 
can be applied directly to the low-x behavior of 
the integrand of the Bjorken sum rule: gi(x,Q 2 ) — 
g™(x,Q 2 ) = \[Au(x,Q 2 )- Ad{x,Q 2 % as well as to the 
other nonsinglet combination Au(x,Q 2 ) + Ad(x,Q 2 ) — 
2As(x,Q 2 ) that also contributes to gf' n (x , Q 2 ) . Note, 
however, that effects similar those in the BFKL pomeron 
may also be important at very low x J35|. The 
flavour-singlet combination of structure functions has 
a more complicated low-cc behaviour, which could be 
important for the extraction of the Aq. More singular 
low- a; behaviours have been proposed in the literature 
motivated by non-perturbative QCD considerations 
|3q, CTl . It is not clear whether the behaviour in 



equation (25) is relevant to the data presently available: 
one SMC data point may be in its region of applicability, 
and could in principle be used to normalize the 
perturbative QCD formula, serving as a basis for 
extrapolating the integrals to x = 0. In practice, it does 
not seem at present that this would have a significant 
effect on the evaluation of the Bjorken sum rule. 

The analysis of the polarized structure function 
data has often assumed that the transverse polarization 
asymmetry 

dal- - da^ (28) 



A i = 



da^~* + da^ 

is negligible. This is related to the spin-flip photon 
absorption asymmetry 

2&TL 



A 2 = ^ (29) 

Cl/2 + °3/2 

and the longitudinal Ai asymmetry (|^) through the 
relation: 



.L = d (A a - 7 (l-|)Ai 



(30) 



were <Ji/ 2 and a^/ 2 are the virtual photon-nucleon 
absorption cross sections for total helicity 1/2 and 
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3/2, respectively, cttl arises from the helicity spin- 
flip amplitude in forward photon-nucleon Compton 
scattering, 7 = 2Mx/-\/Cp, and y = u/Ei epton , 
where v is the energy transfer in the laboratory frame. 
The coefficient d is related to the virtual photon 
depolarization factor D by 



D 



(31) 



The asymmetries A\ and A 2 are subject to the following 
positivity conditions Q 

1^4x1 < 1 , \A 2 \<t/R. (32) 
and are related to the structure functions g± i2 by 



Ai 



1 



(91 



I 2 92) 



A, 



7 



(.91 +.92), (33) 



where Fi = i 7 2(l+7 2 )/2x(l+i?) is the spin-independent 
structure function. Recently, data on g 2 have become 
available for the first time (3^, [h|. They indicate that 
it is considerably smaller than the positivity bound in 
equation (32), and is very close to the leading twist 
formula of ref. lilt]: 



gr(x,Q 2 ) = ~9i(x,Q 2 ) + 



J 9i(t,Q 2 )j, 



(34) 



The data on g 2 are also compatible with the Burkhardt- 
Cottingham sum rule |42| 



g 2 (x, Q 2 )dx — 



(35) 



which has been verified to leading order in perturbative 
QCD. 43, 13. However, the experimental errors are still 
considerable, in particular because the low-x behaviour 
of g 2 is less well understood than that of g\. However, 
the data already tell us that the uncertainty in g 2 is not 
significant for the evaluations of the T\' n . 

1.4-- Higher orders in QCD perturbation theory: 

The perturbation series in QCD is expected to be 
asymptotic with rapidly growing coefficients: 



S(x) 



n=0 



7T 



x = — , c„ ~ n\K n n J 



(36) 



for some coefficients if, 7 |45|, [46|. This type of 
behaviour is associated with the presence of the 
renormalon singularities, as we shall discuss shortly. 
Such series are often evaluated approximately by 
calculating up to the "optimal" order, implicitly defined 
by 

\c nopt x n ^\ < \c nopt+1 x n ^ +1 \ (37) 



and assuming an error of the same order of magnitude 
as c n opt x n ° pt . The question arises whether one 
can approach or even surpass this accuracy without 



calculating all the terms up to order 



^opt • 



This 



possibility has been studied using the effective charge 
(ECH) approach @, g§ and using commensurate scale 
relations |49|, |5(J . In this section, we discuss the use of 
Pade approximants (PA's) for this purpose [pl| p2| . 



Pade approximants 53, p4| are rational functions 



chosen to equal the perturbative series to the order 
calculated: 



[N/M] 



N 



ao + a\X + ... + ajyx 
1 + hx + ... + b M x M 

[N/M] = S + 0{x N+M+1 ) 



(38) 



Under certain circumstances, an expansion of the PA 
in equation ( |3^ ) provides a good estimate, c^_ M+1 , 
the Pade Approximant Prediction (PAP), for the next 
coefficient cjv+m+i in the perturbative series |55| . For 
example, we have demonstrated that if 



c n C n +2 



-n+l 



1 

* 

71 



(39) 



as is the case for any series dominated by a finite number 
of renormalon singularities, then <5[jv/mi defined by 



'[N/M] 



"N+M+l - CN+M+1 

cn+m+i 



(40) 



has the following asymptotic behaviour 



J[N/M] 



Ml 
Tm 



where L = N + M + a M (41) 



and where a is a number of order 1 that depends on the 
series under consideration. This prediction agrees very 
well with the known errors in the PAP's for the QCD 
vacuum polarization D function calculated in the large 
Nf approximation |5^|, as seen in Fig. ||a. 

Large-iVf calculations of the perturbative correc- 
tions to the Bjorken sum rule (5?]] indicate the presence 
of only a finite number of renormalon singularities, so 
that the PAP's should be accurate. Using the known 
terms in equation @, the [1/2] and [2/1] PAP's yield 
the following estimates for the fourth-order coefficient 
l5l: 



-4[PA] 



Bj 

"a[pa\ 



— Ill ([1/2] PA) 



-114 ([2/1] PA) 



(42) 



and the error estimator in equation (41) with a = 1 
yields 

S [1/2] ~ -1/8; S [2/1] ~ -1/4 (43) 
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Figure 2. Relative errors in the [N/M] Pade approximants 
(a) to the QCD vacuum polarization D-funciiDn, evaluated 
to all orders in the largc-Nf approximation |56| (the rate of 
convergence agrees with expectations for a series with a discrete 
set of Borel poles), and (b) to the Borel transform of the 
D-function series, where the convergence is particularly stri king 
The straight lines correspond to the error formulae, eqs. (kll) and 
(fell), respectively. 



These results can be combined to obtain 



C 4[PA] 



-111 

1 + #[1/2] 



-114 

+ £[2/1] 



which is very close to the ECH estimate 



C 4[ECH] 



-130 



-139 (44) 



(45) 



A second application of PA's is to "sum" the full 
perturbative series. The latter is ambiguous if the 
perturbative series possesses an infrared renormalon 
singularity, i.e. a divergence of the form in equation 
(|36|) with K > 0. Consider the following toy example: 



nl x n 

o 



1-xt 



dt 



x .10 



-y/x 



dy (46) 



which exhibits an infrared renormalon pole at y = 1. 
One possible way to define the ambiguous integral on 
the right hand side of equation ( [l6"| ) is via the Cauchy 
principle value prescription [pSf . We see in Fig. | that 
the errors in the Pade "Sums" (PS's) [N/M](x) are 
smaller than the truncated perturbative series 



N+M 



(47) 



when n < n opt , which is 5 in this example. You will 
notice in Fig. || that the errors in the PS's become 



unstable for large n: this is because of nearby poles 
in the denominator of equation ( p8| ) which are not 
important for small n. Also shown in Fig. || as 
"combined method" is a systematic approach to treating 
these poles and optimizing the PS's for large n, which 
is described elsewhere fr9l. 
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Figure 3. The relative errors between partial sums of the 
series S(x) = Sn!i™ and the Cauchy principal value of the series 
(solid line) is compared with the relative errors of Pade Sums 
(dotted line). We see that the relative errors of the Pade 
Sums are smaller than those of the partial sums in low orders, 
fluctuate in an intermediate regime, and are again more 
accurate than the partial sums in higher orders. The 
fluctuations are associated with nearby poles in the Pade Sums, 
that may be treated by the "combined method" mentioned in 
the text, shown as the dashed line. 



3— rd order partial sum 

4— th order partial sum (ECH) 
- - [1/2] Fade 
[2/2] Pade 




1.25 



Figure 4. The scale dependence of cr s (3GeV 2 ) obtained from a 
fixed value f(x) = (6/ff^) X Q.164 = 0.783, (cf. eq. (g) ) , for 
Q/2 < /i < 2Q, using the naive third- and fourth-order per- 
turbative series and the [1/2] and [2/2] PS's. 



Evidence that the PS's for the Bjorken sum rule 
provide a good estimate of the perturbative correction 
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factor in equation ( |iq ) is provided by the study of the 
renormalization scale dependence. We see in Fig. || 
that the renormalization scale dependence of the [2/2] 
PS is much smaller than that of the [2/1] and [1/2] 
PS's, which is in turn much smaller than that of the 
naive perturbation series evaluated to third order. We 
recall |5(| that the full correction factor should be scale- 
independent, and interpret Fig. [|as indicating that the 
PS's may be very close to the true result. 



in this plane 



f(x) 




[3/3] Pade - ECH — 



0.000 0.025 0.050 0.075 0.100 0.125 
x=a V 7T 



Figure 5. Different approximations to the Bjorken sum 
ruie correction factor f(x), third-order and fourth-order 
perturbation theory, [1/2], [2/1] and [2/2] Pade Sums are 
compared. Also shown as a vertical error bar is the val ue o f f(x) 
we extract from the available polarized structure data (p4|). 

Fig. H shows the estimates of the perturbative 
QCD correction to the Bjorken sum rule obtained in 
various different approximations, including third-order 
perturbation theory, fourth-order perturbation theory 
estimated using the ECH technique and the [2/1], [1/2] 
and [2/2] PS's. We interpret the latter as the best 
estimator, and take the difference between it and the 
[2/1] and [1/2] PS's as a theoretical uncertainty. Also 
shown in Fig. |B| is the experimental error on this 
quantity, as extracted from the combined analysis of 
the available experimental data discussed in the next 
section. 

More information can be extracted by considering 
PA's in the Borel plane. The Borel transform of a 
perturbative series is defined by 



S(x) = c n x n 



Borol 



C n +1 



n+1 



(48) 



where (3o = (33 — 27V/)/3. A discrete set of renormalon 
singularities would show up as a set of finite-order poles 



{y-ykY 



(49) 



The PA's in equation (j^) are clearly well suited to find 
the locations yu and the residues of such poles. In the 
case of a perturbative series dominated by a finite set of 
L renormalon singularities, a sufficiently high-order PA 
will be exact 



[M/N] (y) = S(y) : for M + N > L 



(50) 



for some Lq oc L. Generically, in any case where the 
quantity analogous to eq. (p9|), e„ ~ l/« 2 , the error 
analogous to (40) is given by 



n 2 



>[M/M] 



(Ml) 

L 2M 



(51) 



This prediction of very rapid convergence is confirmed in 
Fig. ||b |5l[] in the case of the QCD vacuum polarization 
D function evaluated in the large Nf limit (56), which 
has an infinite number of renormalon poles. 



y- [2/1] Fade, with pole at y - 1.05 
fit with poles at y = -1. 1. 2 
□ fit with poles at. y = -1, 1, 2, 3 
-t tit with poles at y = -1, -2. 1, 2 
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Figure 6. The locations and residues of poles in the [2/1] PA 
and in rational-function fits to the Borel transform of the first 
four terms in the perturbation series for the Bjorken sum rule. 
We note that the location of the lowest-lying infrared renormalon 
pole is estimated accurately by Pade Approximants in the Borel 
plane, and that its residue is stable in the different fits. 

The power of the application of PA's in the Borel 
plane is shown in Fig. [g, where we see that the [2 / 1] PA 
of the Borel-transformed Bjorken series has a pole at 
y = 1.05 . The agreement with the expectation of a first 
infrared renormalon pole at y = 1 is striking. Fig. |^ 
also shows other rational-function approximations to 
the Borel-transformed Bjorken series. We see evidence 
that this is dominated by a strong infrared renormalon 
pole at y = 1, that there may be a weak ultraviolet 
renormalon pole at y = — 1 and possibly another pole at 
y = 2. 
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The ambiguity in the definition of the perturbative 
Bjorken series associated with the y — 1 renormalon 
singularity corresponds to a possible l/Q 2 correction of 
magnitude |52] 



A(rf-r?) -±^o.98^^ 

6 Q z 



(52) 



It is expected that the QCD correction to the Bjorken 
sum rule should include a higher-twist correction of 
similar form with magnitude pi] , [62] , |63[ | 



Table I 



pP,"! curren tly available from experiments at CERN and SLAC. 



experiment 


target 


Ti 


E142 


n 


-0.045 ±0.009 


E143 


P 


0.124 ±0.011 


E143 


d 


0.041 ± 0.005 


SMC 


d 


0.023 ±0.025 


SMC ('94) 


d 


0.030 ±0.011 


SMC 


P 


0.122 ±0.016 


EMC 


P 


0.112 ±0.018 



iHT(r; _ r;) _o™ (a) 



The perturbative ambiguity in equation (|52j) is cancelled 
by a corresponding ambiguity in the definition of the 
higher-twist contribution. In the next section we will 
treat equation (|5^) as a correction (with error) to be 
applied to the perturbative QCD factor shown in Fig. ||. 

We have also compared PA's to the predictions of 
commensurate scale relations within the framework of 
ref. ^0|. The predictions of the two approaches are 
numerically very similar, and we give formal reasons in 
ref. [|59] why we believe that this should be so. However, 
we shall not use commensurate scale relations in the 
data analysis of the next section. 



1.5. Numerical analysis of the Bjorken sum rule 

Table I shows the data on the integrals Pf'™ currently 
available from experiments at CERN and SLAC 
|p7|,||64ft-[pq|. We do not attempt to correct these 
numbers for any of the effects discussed in section 1.4, 
such as the Q 2 -dependence of the asymmetry Ai, the 
extrapolation to low x, or the transverse polarization 
asymmetry. We do not believe that any of these effects 
will change any of the data outside their quoted errors. 
We choose to evaluate the Bjorken sum rule at Q 2 = 3 
GeV 2 , which requires rescaling all the data as described 
in ref. il. 



All experimental data have been evolved to Q 2 = 3 GeV 2 . 

The following is the combined result that we find for 
the Bjorken sum rule: 

rP(3GeV 2 ) - r^SGeV 2 ) = 0.164 ± 0.011 (54) 

which is indicated by a vertical error bar in the lower 
left corner of Fig. |[ Comparing this value with the [2/2] 
PS estimate also shown there, we find 



a s (3 GeV 2 ) = 0.328 



+0.026 
0.037 



which becomes 
a s (M 



0.119 



+0.003 
0.005 



± 



(55) 



(56) 



when we run a s up to Af| using the three-loop 
renormalization group equation. The errors quoted 
in equations (|5^) and ( p6[ ) are purely experimental, 
and the second ± sign in equation (^6|) indicates that 
further theoretical systematic errors must be estimated. 
Those we have evaluated include that associated with 
the renormalization scale dependence shown in Fig. ^ 
(±0.002), the difference between the [2/2] and [2/1], 
[1/2] PS's (±0.002), and the correction due to the 
higher-twist estimate in equation (|53| ) (—0.003 ± 0.002), 
whereas the error in the running of a s is found to be 
negligible. Combining these estimates with equation 
©, we find H 



a s (M 2 z ) = 0.116±°;°<g ± 0.003 



(57) 



The stability of this result is indicated in Fig. ^, 
where we exhibit the values of a s (Af|) obtained using 
different orders of perturbation theory, compared with 
our result obtained using the [2/2] PS. Also 

indicated is the shift induced by the higher-twist 
correction, which lies within our error bars. 



S 



1.6. Decomposition of the Nucleon Spin 



8™ 



0.14 



0.13 



0.13 



X usual P.T 

[1/2] & [2/1] Pade 
+ [2/2] Pade 

1 2 ; 



Figure 7. Values of a s (M^) obtained using different orders of 
perturbation theory, compared with our result (p6|), obtained 
using the [2/ 2] The size of the shift induced by higher-twist 

correction is ( |53|) indicated by a downward arrow to the right of 
the [2/2] point. 



average: 
0.27 ± 0.04 



X = 2." 



-0.2 

parton 0(o, ) 0(a^ ) 0(t^ ) 0(dj ) 

model 

E142 + E143-p • E143-d <> SMC-d(92) S SMC-d(94) □ SMC-p S EMC 
all A£ values evolved to Q 2 = 3 GcV s 



BjSR 
GLSB 
R r 

pp. pp -* 7+i cts 

LGT(T-system) 

ce, bb decays 

scaling violation (DIS) 

jet. rates (HERA) 

pp -* bb+jets 

H 7 

scaling violation (e + e ) 
event shapes (e e ) 
pp -> W+jets 
LEP, R 7 



p (GoV) 
1.73 
1.73 
1.777 
4.0 
5.0 
6.6 
7.1 

18.0 

20.0 

33.0 

50.0 

60.0 

80.2 

91.2 



0.075 0.100 
« S (M 7 ) 



Figure 8. Compilation of world . 
different sources (adapted from Ref. 1 



ita. 
1). 



from 



As can by seen from the compilation in Fig. g, our 
central value for a s (M§) is quite compatible with other 
determinations and the world average, which is quoted 
to be § |7|, !§ 



<Xa{M%) 



0.117 ±0.005 



(58) 



Indeed, the error quoted in equation (^7J) is quite 
competitive with the most precise determinations of 
a s (M§) that are available. Moreover, we note that 
plenty of precise new data will soon be available from 
the SMC experiment at CERN, the E154 and E155 
experiments at SLAC, and the HERMES experiment 
p9t at DESY. In the longer run, experiments with a 
polarized proton beam at HERA will provide valuable 
information on the behaviour of g\ at low well as 
on its CJ 2 -dependence at fixed x. 



Figure 9. The values of A£(Q 2 =3 GeV 2 ) extracted from each 
experiment, plotted as functions of the increasing order of QCD 
pertu rba tion theory used in obtaining AS from the data (from 
Ref. [B2J, updated with most recent data). 

So far, we have only discussed the combination 
— r™ which enters in the Bjorken sum rule. The 
individual T\ ,n can be used as in equation (H), though 
not forgetting the perturbative QCD corrections in 
equation fll7|), to extract the individual Aq. As is 
seen in FigTM the different experiments on both proton 
and neutron targets are all highly consistent, once the 
perturbative QCD corrections are taken into account. 
Some time ago, it appeared as if the neutron data 
from E142 might be at variance with the other data 
points. However, this is no longer the case if all 
the higher-order corrections in equation ( |l7| ) are taken 
into account, and the latest evaluations of the 
E142 data indicate a different preliminary value of 
T", as seen in Table I. Making a global fit, we find 



Au = 0.82 ± 0.03 ± 
Ad = -0.44 ± 0.03 ± 
As = -0.11 ± 0.03 ± 



and 



AE = Am + Ad + As = 0.27 ± 0.04 ± . . 



(59) 



(60) 



where the second ± sign indicates that further 
theoretical and systematic errors remain to be assigned. 
These include higher-twist effects, errors in the 
extrapolation to low x which is more complicated than 
for the nonsinglct combination of structure functions 
appearing in the Bjorken integrand, the possible Q 2 - 
dependence of A 1: etc.. We believe that these errors 
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may combine to be comparable with the errors quoted 
in equations ( |59| ) , ( |60| ) , but prefer not to quote definitive 
ranges for the Aq until all these errors are controlled as 
well as those appearing in the Bjorken sum rule. 

A£ vs As 



0.4 - 



i 




x E142 4- E143-p o E143-d SMC-d(92) s:s SMC-d(94) □ SMC-p X EMC 

Figure 10. The values of AE and As extracted from each 
experiment, plotted against each other. All data have been 
evolved to common Q 2 = 3 GeV 2 . The clear linear correlation 
between AS and As results from the linear relations 

One may also get a feeling for the expected range 
of AS and As by plotting the results for these 
two observables extracted from each of the existing 
experiments, as shown in Fig. [l^. 

2. Outlook 

In this talk we have concentrated on the phenomenolog- 
ical analysis of the data on polarized structure functions 
presently available. As we have seen, these tell a remark- 
ably consistent story, once higher-order QCD correc- 
tions are included. We have not addressed in great detail 
here the theoretical interpretation of the data, nor their 
spin-offs in hadron physics and elsewhere, nor possible 
future developments in this field. In fact, these mea- 
surements provide valuable insights into important is- 
sues in non-perturbative QCD, such as the role of chiral 
symmetry in nucleon structure [|Ts(] , the axial anomaly 
and the £7(1) problem [^o)-|2^], and the relationship be- 
tween current and constituent quarks [f70f-|f75| which are 
provoking lively theoretical debates (see Ref. Q for a 
recent application to the pion structure). 

The polarized structure function data support previ- 
ous indications from the 7r-nucleon cr-term and elsewhere 
that strange quarks in the nucleon wave function cannot 
be neglected, with interesting implications for the analy- 
sis of recent data from LEAR on <fi production in proton- 
antiproton annihilation |77| . Among other spin-offs, we 
recall that the axial-current matrix elements extracted 



from polarization data determine scattering matrix el- 
ements for candidate dark matter particles such as the 
lightest supersymmetric particle |7^] and the axion |32j . 

In the future, we look forward to the completion of 
the SMC programme and its possible HMC successor at 
CERN, the E154 and E155 experiments at SLAC, data 
from the HERMES experiment at HERA, the polarized 
proton programme at RHIC, and possible polarized 
electrons and protons in the HERA ring. The tasks of 
these experiments will include the determination of the 
gluon contribution to the nucleon spin, the elucidation 
of the Q, 2 dependence of A\, and the low- a; behaviour of 
<7i . These will continue to fuel activity in this interesting 
field for the foreseeable future, which will lead us to 
a deeper understanding of the nucleon, an object we 
thought we knew so well, but which reveals a new face 
when it spins. 
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